Dynamical shift conditions for the Z4 and BSSN formahsms 
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A class of dynamical shift conditions is shown to lead to a pseudo-hyperbolic evolution system, 
both in the Z4 and in the BSSN Numerical Relativity formalisms. This is done by using a plane- 
wave analysis which can be viewed as an extension of the standard Fourier analysis for this kind 
of systems. The proposed class generalizes the harmonic shift condition, where light speed is the 
only non-trivial characteristic speed, and it is contained into the multi-parameter family of minimal 
distortion shift conditions recently proposed by Lindblom and Scheel. The relationship with the 
analogous 'dynamical freezing' shift conditions used in black hole simulations discussed. 

PACS numbers: 



I. INTRODUCTION 

General covariance is a characteristic property of Ein- 
stein's theory of Gravitation. There are four coordinate 
degrees of freedom in the field equations, allowing to 
freely choose the spacetime coordinates or, in the 
framework of the 3-f 1 decomposition, the lapse function 
a and the shift vector /3'. 

This gauge freedom can be used, like in the early years 
of General Relativity, to set up a complete evolution sys- 
tem (consisting of the field equations plus the gauge con- 
ditions) with a well posed Cauchy problem. The well 
known harmonic coordinate conditions 



0, 



(1) 



provide a simple way to get a well posed initial value 
problem P| with an extremely simple principal part, con- 
sisting of one wave equation for every component of the 
four-dimensional metric tensor: 



(2) 



0, 13 , where the dots stand for terms not belonging to 
the principal part. 

Although the resulting system is still currently used 
in analytical approximations, its use in Numerical Rela- 
tivity is very limited, mainly because the four conditions 
completely exhaust the gauge degrees of freedom, and 
there is no flexibility left that could be used to fit the pe- 
culiarities of the specific systems one wants to model (but 
see also some generalizations in Refs. Q to 0)- 

The current alternative, represented by the 'new hy- 
perbolic formalisms' (Refs. to 0|), is to use somehow 
the momentum constraint as a tool for ensuring hyper- 
bolicity, instead of the three space coordinate conditions 
in 10. In this way, only the time gauge condition 



Dx° = 



(3) 



is kept, or one of its generalizations (harmonic slicing). 
This happens to be very convenient for numerical simu- 
lations, because © implies a direct relationship between 



the lapse a and the spatial volume element ^/7, which can 
be used to avoid collapse singularities . The main ad- 
vantage is that one can use now the shift degrees of free- 
dom to simplify the system, either using normal coordi- 
nates (zero shift) or any other kinematical choice adapted 
to the specific problem under study. 

Recently, there has been a renewed interest in the use 
of dynamical shift vectors. In the context of the BSSN 
formalism j23, 0| , some shift conditions have been pro- 
posed [i^l that manage to 'freeze' black hole dynamics 
near the horizon, leading to long term numerical simula- 
tions although. As we will see later, the harmonic shift 
condition given by is similar to (but not contained in) 
the ones discussed in Ref. ^2^. 

It is clear that the principal part of the original BSSN 
system is modified by the choice of a dynamical shift, al- 
though no hyperbolicity analysis of the modified system 
has been yet published, to the best of our knowledge. 
This is not surprising because the BSSN formalism, like 
the ADM one, is of a mixed type: first order in time, 
but second order in space, and therefore the standard 
Fourier analysis would lead to the conclusion that the 
mixed order system is not hyperbolic (2^] . This has been 
explicitly shown by Fritelli and Gomez for both the ADM 
and BSSN formalisms We will present here an al- 

ternative plane-wave analysis, based on the underlying 
physics, in order to reveal a related property, which we 
will call 'pseudo-hyperbolicity' to avoid confusion. As far 
as the underlying physics does not change when passing 
from the fully second order system to the mixed order 
version of the same equations, pseudo-hyperbolicity can 
be seen as the imprint left on the mixed order system by 
the true hyperbolicity of the fully second order version 
which was at the starting point. 

From a different point of view, a generalization of the 
minimal distortion shift condition has been proposed by 
Lindblom and Scheel in the context of the first order 
KST formalism p^ . Surprisingly enough, in spite of the 
fact that the resulting system contains at least twenty- 
two free parameters, none of the cases discussed in |25j 
verifies the condition that all the non-trivial character- 
istic speeds do coincide with light speed. The surprise 
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comes from the fact that one would expect this condition 
to be ensured by the use of the full set of harmonic co- 
ordinates (0, which are actually a particular case of the 
gauge conditions proposed in 25]. 

As a contribution to clarify this issues, we will present 
here a family of dynamical shift conditions, which is 
closely related with the ones presented in '25| and [25l| . 
We will do so first in the framework of the general- 
covariant Z4 formalism |^ , extending then the results to 
the BSSN c ase, which can be derived from the Z4 one in a 
simple way \2T\ . We will perform a complete plane- wave 
analysis of both systems. The resulting characteristic 
speeds are directly related with the main free param- 
eters of the proposed family. The further requirement 
that all the non-trivial characteristic speeds do coincide 
with light speed, allows one to recover the harmonic case. 



II. THE Z4 SYSTEM 

A. The Evolution Equations 

The Z4 covariant formalism introduces a four- 
dimensional vector as a supplementary dynamical field 
Zfj,. The evolution equations are obtained by adding the 
(symmetrized) covariant derivatives of to Einstein's 
field Equations: 

i?p. + V^Z, + V,Z^ = 8 7r(T^,-iTg^,) . (4) 

In the 3-1-1 decomposition the line element is written as: 



ds^ = -a^ df + 7,y {dx' + dt) [dx' + (3' dt) (5) 

where a and /?* are the lapse and the shift, respectively, 
and 7ij is the spatial 3-metric. Using this decomposition, 
the general covariant equations Q do consist of a system 
of pure evolution equations: 



{dt - Cp) jij -2 a Kij 

- 2Kf. + {trK - 20) K,j 



(6) 



- ^,, + -(tr^- r)7.,] (7) 
{dt -Cp)e = I f^R + 2 VkZ^ + {ti-K - 2 9) tiK 

- tiiK^) -2 Z''ak/a-2T] (8) 
{dt - Lp) Z, = a [V, {K^ - 5,HtK) + d,Q 

- 2K^ Z,-Qa,/a-S^] (9) 

where we have noted 

e = a Z°, r = STra^ T™, S, = 87ra r°, S.,j = Stt . 

(10) 



In a recent work |27|, a symmetry breaking mecha- 
nism is proposed that, starting from the Z4 system H6I9|I . 
allows one to recover an evolution system which is equiv- 
alent, up to quadratic source terms, to the BSSN system 
12(1 l2ll| (partial symmetry breaking). Also, in the first 
order case, the same mechanism allows one to recover the 
multi-parameter KST system [1Z| (full symmetry break- 
ing) orjto be more precise, a 'live gauge' version of the 
same [l8|. 



B. Gauge evolution equations 

The harmonic gauge conditions can be easily ex- 
pressed in the 3-1-1 formalism as 



{dt - P^'dr) a = - o? \xK 



(11) 



{dt - d'-dr) 



[dM(^^)+ d,f^], (12) 



where the first equation is the (harmonic) slicing con- 
dition, whereas the second one provides the (harmonic) 
shift once the slicing is known. 

The harmonic slicing condition Hll|) has been general- 
ized in the context of the Z4 system as follows I27t|: 

{dt- P^'dr) a = - f [tiK -mO], (13) 

where the parameter / is directly related with the gauge 
propagation speed, whereas m provides a coupling with 
the energy-constraint-violating modes, represented by 
the quantity 9. We will see in the following analysis that, 
if one wants the gauge speed to coincide with light speed 
(/ — 1), then a pseudo-hyperbolic system is obtained 
only if m = 2, so the coupling given by the m parameter 
can not be neglected. This conclusion coincides with the 
result of Ref. 27], where it was confirmed by the robust 
stability numerical test. 

The harmonic shift evolution equation (|12|l can be gen- 
eralized along the same lines: 

{dt-p^'dr) = ~a^ [2fiV' + adHna-ddHny/^~'nl3' 

(14) 

where we have defined 



(15) 



Notice that the advection term on the left-hand-side, 
which was absent in Ref. is needed if one wants 

to recover the harmonic shift as a particular case. As 
we will see in the following analysis, the parameters /i 
and d are directly related with the characteristic speeds 
of the longitudinal and transverse shift components, re- 
spectively, in the same way as / is related with gauge 
speed. The parameter a, instead, has no direct relation- 
ship with the characteristic speeds: its role is very simi- 
lar to the parameter m in the lapse condition, as we will 
see that specific values of a will be required to ensure 
pseudo-hyperbolicity in degenerate cases, so one can not 
just neglect this kind of coupling. The parameter 77, in 
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turn, corresponds instead to a damping term which has 
shown to be crucial to get stable long term simulations 
|22| |. We have not included, however, the rj term in our 
analysis to avoid masking the genuine wave propagation 
effects with the artificial damping produced by this kind 
of terms. 



III. LINEAR PLANE- WAVE ANALYSIS 

The system ()6I9|I is of a mixed type: first order in 
time but second order in space. This means that, accord- 
ing to the standard methods j2^, based on the Fourier 
analysis of the principal part, it can not be classified as 
hyperbolic. This is also the case of the original ADM 
system |23|, where the quantities (6, Zk) are supposed 
to be zero. In what follows, we will present an alter- 
native plane-wave analysis, starting with the ADM case 
first and including the supplementary quantities (0, Zk) 
later. 

It is well known that any metric can be written down 
at a given spacetime point P in a locally inertial coordi- 
nate system such that the first derivatives of the metric 
coefficients vanish at P. We will take advantage of this 
to write down the line element at P as 



-al dt^ + 7°,- {dx' + (3f} dt) {dx^ 



(16) 



It is clear that the validity of the expression 1)16(1 is 
strictly local: second and higher order derivatives of the 
metric coefficients at P can not be supposed to vanish: 
they are rather related one another by the field equa- 
tions. This suggests the splitting of the metric into two 
components: 

• A uniform static background of the form 1)16(1 

• A dynamical perturbation which, when superim- 
posed to the background in a linear way, allows 
one to recover the full metric. 

It makes sense then to decompose the dynamical per- 
turbation into plane waves, with a space dependence 
given by 



ao ~ e 



UJ-X ~ / 



(3^ 

1^3 - 



a{uj, t) 



(17) 
(18) 
(19) 



where lou — ujuk, — 1. 

Up to here, we have followed the standard Fourier anal- 
ysis. Now we will depart from the standard path by de- 
composing the dynamical variable Kij in a form which is 
consistent with the exact evolution equation ©, namely 



(20) 



where one must notice the i lo factor on the right-hand- 
side. This extra factor is not present in the standard hy- 
perbolicity analysis, where only the principal part of the 



system is used, breaking in that way the direct relation- 
ship (|SJ) between the metric and the extrinsic curvature. 

Note that equation ((HJ is crucial to relate the original 
(second order in time) version of the field equations Q 
with the resulting (first order in time) 3-1-1 version ((619(1 . 
This is why we will choose the alternative decomposition 
((2(JII in order to look for the imprint in the 3-1-1 system of 
the hyperbolicity properties of the original second order 
version. 

The same thing can be done with the supplementary 
quantities, which can be considered as an additional per- 
turbation of the Einstein's equations background. The 
original system Q is of first order in Z^, but it can be 
viewed alternatively as being of second order in some 
'potential' quantities which time derivative can be de- 
fined to be precisely Z^ . In this way l@J could be seen as 
a fully second order system in (7^,^, Y^) and the same ar- 
guments as before would justify the following plane- wave 
decomposition of the supplementary quantities 

e = iuje''^-''Q{uj,t) (21) 
Zk = lue'^-^ Zk{uj,t), (22) 

where the i lo factors still appear on the right-hand-side. 



IV. PSEUDO-HYPERBOLICITY OF THE Z4 
SYSTEM WITH DYNAMICAL SHIFT 

We will perform here a linear plane-wave analysis of 
the Z4 system. This means to substitute the perturba- 
tions described in the preceding section into the evolution 
equations ((6191 II 31 IT^ . keeping only the linear terms. We 
get: 

iw/^o") a 



{dt — i,uj 



-iualfitrK-mQ] (23) 
—i uj [(2 PL ~ d)n} tr7 
2^(7™ + Z')+a7i*(5/ao] (24) 

(25) 

[dt-ioj /?(,") 8 = -iujao [tr7 - 7"" - Z"] (26) 
{dt-iLo /3o") Z, ^ -iLo ao [n, (tiK - 6) - K^] (27) 
[dt - iuj /3q") k^j = -iuao [j^j + Ui nj (tr7 -f a/ao) 

- fli {jnj + Zj) - Uj {%i + Zi)] (28) 

where the letter n replacing one index means contracting 
that index with n^. 

Notice that we have kept the linear source terms in 
((231 125() , in contrast with the usual practice in the stan- 
dard Fourier analysis, where only the principal part of 
the system is considered. In fact, our plane wave analy- 
sis includes (up to the linear order) all the source terms 
(with the only exception of the artificial damping one for 
the shift, as discussed before). Therefore, the underlying 
physics is accounted for in a consistent way. In partic- 
ular, the direct relationship between the metric and the 
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extrinsic curvature is fully preserved. This means that 
the characteristic speeds we are going to compute should 
be the same ones that could be obtained from either the 
fully second order or the fully first order versions, where 
the standard Fourier analysis can be applied in a way 
which is consistent with the underlying physics of the 
problem. 

The system can be written in a compact way 

as 

dtu = -iLu[A- l3fP I]u, (29) 

where u is the perturbation array. The geometric prop- 
erties of matrix on the right-hand-side (characteristic 
matrix) are obviously related with the dynamics of the 
plane-wave perturbations. Allowing for the trivial struc- 
ture of the shift term, the pseudo-hyperbolicity of the 
evolution system (|29|) will depend on the properties of 
the main matrix A. We will say that the system (|29|) 
is 'pseudo- hyperbolic' if and only if A has real eigen- 
values and a complete set of eigenvectors, that is, the 
number of independent eigenvectors must be the same as 
the number of independent dynamical fields (20 in our 
case). This is the analogous of the strong hyperbolic- 
ity property of first order systems. The use of the term 
'pseudo-hyperbolicity' is just to avoid confusion. 

We can start computing the eigenmodes which do not 
contain shift terms. In this list we have: 

Energy cone There are 2 9-related eigenmodes, prop- 
agating with light speed, v = — /Sg" i ao : 

e ± {tn - 7"" - ^") • (30) 

Light cones There are 10 more eigenmodes propagating 
with light speed, v = — /3o" i ao : 

(31) 

Kab ± lab, (32) 

where the letters a, b replacing an index mean the 
projection orthogonal to Ui. 

Lapse cone There are 2 a-related eigenmodes propagat- 
ing with speed v — —(3q ± ag \/J : 

V7 [trK - Bi 9] ± [a/ao + {2 - B,) {tn - 7™ " ^")] , 

(33) 

where have used the shortcut 

mf-2 

Bi ^ -j^ ■ (34) 

The factor / must be greater than zero for pseudo- 
hyperbolicity. Notice that, in the degenerate case 
/ = 1 (harmonic slicing), a well defined pair of 
eigenmodes is obtained only if to = 2, so that the 
parameter B\ can take any value (arbitrary mixing 
with the energy cone). 

The shift-related cones are: 



Transverse shift cones There are 4 eigenmodes prop- 
agating with speeds v ~ — /3o" i Q^o \/m i 

0a/ao)±2^{jna + Za). (35) 

The factor fi must be greater than zero for pseudo- 
hyperbolicity. Notice that, in the vanishing shift 
case, they reduce to the second term, which would 
correspond to standing eigenmodes (zero character- 
istic speed). 

Longitudinal shift mode There are 2 eigenmodes 
propagating with speed v = —Pq ± ao Vd, 

Vd [ p"/ao + B2trK + B3e]± 
[ {a + B2) a/ao — dtrj 
+ (2Ai + 2B2 + B3)(tr7-7""-Z")] (36) 

We have used again the shortcuts 

_d-af _ {2-mf)B2 + 2fi-amf 
^2 = 73^-' = . (37) 

A necessary condition for pseudo-hyperbolicity is 
that the factor d should be greater than zero. This 
condition is also sufhcient in the generic case where 
d is different from both 1 and /. Notice that in the 
degenerate cases one would need to impose addi- 
tional conditions on the free parameters in order to 
get a well defined pair of eigenmodes. For instance, 
ifd — f one must have a = 1, so that the parameter 
B2 can take any value (arbitrary mixing with the 
gauge cone) . If we have further degeneracy, that is 
d = f = 1 (remember that / = 1 implies to = 2), 
then it follows from H37|) that /i = 1 also, so that 
one gets the harmonic shift case. 

In summary, there are 20 fields in the evolution sys- 
tem and we have got real characteristic speeds and 20 
independent eigenvectors, provided that all the charac- 
teristic speed parameters /, /i, d are greater than zero. 
The system in then pseudo-hyperbolic in the generic case, 
although degenerate cases, where different characteristic 
speeds actually coincide, require additional conditions on 
the remaining parameters a, to. 

V. PSEUDO-HYPERBOLICITY ANALYSIS FOR 
THE BSSN SYSTEM WITH DYNAMICAL SHIFT 

A pseudo-S£ectral analysis of the original BSSN 
system |2^ has been done recently p9| . We will 
proceed here instead with the linear plane-wave analysis 
of the complete system, including the dynamical shift 
terms. We can take advantage of the symmetry breaking 
mechanism proposed in Ref. ^7}. Starting from the Z4 
system equations (|nito|3l, we will take the following steps: 

1. Perform the dynamical fields recombination 

TL 

K\,^K,,--ei,, (38) 
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2. Suppress 9 as a dynamical quantity, setting its 
value equal to zero wherever it appears in the evo- 
lution equations. 

This process alters the evolution equation for the ex- 
trinsic curvature Kij, even if one has Kij = K'ij after 
the second step. One gets as a result a one parame- 
ter family of evolution systems, with different principal 
parts for every value of the n parameter, namely: 

{dt - P^'dr) a = - ftxK (39) 

{dt - P"" dr) 13' = - [2^iV' + adUna-^ dd'lny/^ 

- (40) 

{dt - P'^dr) 7y ^ -2 a 

+ l^k{^,p'')+-f,t{^^p'') (41) 

{dt - 13'- dr) + dk{a {8\ trK ~ K^i)] = ■■■ (42) 

{dt - P'^dr) -I- dk[a A^,] = • • • (43) 

where we have noted for short 

2 A^, = + S'l {d.j Ina + dj In^ + 2 Vj) 

+5^^ (9, Ina + d, In^ +2V,)~nV^ 7^ (44) 

In particular, it has been shown in |22| that the prin- 
cipal part of the system obtained when n — 4/3 can be 
rewritten, by further rearranging the dynamical fields, as 
that of the BSSN system. We will name this system Z3- 
BSSN to avoid confusion. Both systems are then linearly 
equivalent (equivalent up to quadratic source terms, see 
Ref. (3(|), so that showing pseudo-hyperbolicity for the 
Z3-BSSN system, as we will do in the present section, 
amounts to show the same for the original BSSN system. 

As far as the second step of the symmetry breaking 
procedure suppresses the dynamical field 9, the linear 
plane-wave analysis must be repeated from scratch, al- 
though the results of the preceding section provide a use- 
ful guide. We obtain the following list of eigenvectors and 
eigenvalues: 

Standing mode There is 1 mode propagating along the 
normal lines, that is with speed v = —I3q' : 

tr^ „ _ i« . (45) 

This is what is left of the Energy cone in the Z4 
system after suppressing 9 as a dynamical field. 

Light cones There are 10 modes propagating with 
speed V = — /3o" i ao, namely 

Kab ± [lab - \ {tn - 7™ - fab] (46) 

Kna ± (47) 
Lapse cone There are 2 modes propagating with speed 

^trK ± a/ ao (48) 



Transverse shift cones There are 4 modes propagat- 
ing with speed v = —(3^ ± ag i 

{(3a/ao)±2^{%,a + Za) (49) 

Longitudinal shift cone There are 2 modes propagat- 
ing with speed v = — /3g ± ao ^fd'- 

Vd[i3''/ao + B2 trK]± 

[{a + B2) 5/ao - dtn + 2 (<r7 - 7"" - ^")](50) 

where the parameter B2 is the same that appears in the 
Z4 system, as defined in (|37|) . Notice that, in the de- 
generate case d = f, a, well defined pair of eigenmodes 
is obtained only if a = 1, so that B2 can take any value 
(arbitrary mixing with the lapse cone). 

In summary, there are 19 fields in the evolution system 
and we have got real characteristic speeds and 19 inde- 
pendent eigenvectors, provided that all the characteris- 
tic speed parameters /, fi, d are greater than zero. The 
system is then pseudo-hyperbolic in the generic case, al- 
though the degenerate case d — f requires the additional 
condition a = 1. Notice that the harmonic case is recov- 
ered precisely when 

= / = /i = 1 . (51) 

Comparing with Ref. 22], the only relevant term in 
their "Gamma driver" shift conditions is the one. It 
is clear that this corresponds to our parameter choices 

a = 0, M=^rf (52) 

so that pseudo-hyperbolicity is ensured provided that 
d ^ f . The main difference, as stated before, is that 
our shift conditions (|14|l are a generalization of the har- 
monic ones (|12ll . This is only possible by keeping, as 
we do, the advection term in the shift evolution equa- 
tion H14|) . This term was suppressed in Ref. 22] in order 
to freeze the dynamics in the neighborhood of the black 
hole's apparent horizon. 

VI. CONCLUSIONS 

In this paper, we have studied a multi-parameter fam- 
ily of dynamical gauge conditions ((T^ ITU) , which gener- 
alizes the harmonic gauge conditions Q along the ways 
sketched in Refs. |23], |2a] (Gamma-driver shift condi- 
tions). Starting with (the second order version of) the 
general covariant Z4 formalism, we have computed ex- 
plicitly all the eigenmodes, identifying the choices of the 
gauge parameters {/, to, /i, a, d} that make the full evolu- 
tion system pseudo-hyperbolic. The relationship between 
the gauge parameters and the characteristic speeds is di- 
rect and simple. Depending on the parameters choice, 
gauge (lapse and shift) propagation speeds can be made 
to coincide or not with light speed. 
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The same kind of analysis has been done with the cor- 
responding gauge conditions, for the second order system 
(EHl-ESJiWhich is Hnearly equivalent to that of the BSSN 
system "30]. The Gamma-driver shift condition which 
has been used in Ref. corresponds to the parameters 
choice H52(l . but with the advection term in the left-hand- 
side of the shift equation H4U|I suppressed. One could redo 
the analysis without that term, just by adding it to the 
right-hand-side of the same equation. This would affect 
the shift cones (|49l I50|l . which would then intersect the 
other ones, leading to a more involved causal structure. 
We have chosen to keep instead the shift advection term 
in place, so that the original harmonic gauge condition 
(|T|) is kept inside our family. 

Finally, we will briefly discuss the relevance of our re- 
sults in connection with the ones presented in Ref. l25j, 
in the context of the (first order) KST formalism. In 
that paper, a 22-parameter family of gauge conditions 
was presented which contains \^ as a sub-family, 
but no parameter choice was found ensuring both sym- 
metric hyperbolicity and the additional requirement that 
all the non-trivial characteristic speeds should coincide 
with light speed. Is this because these two conditions are 
actually incompatible in the KST formalism... or it is 
rather because there are only few 'good' choices (maybe 
just one), hidden in the huge 22-parameter space, just 
waiting to be identified?. 

We can not fully answer this question here because we 
are dealing in this paper just with pseudo-hyperbolicity, 
not with the stronger condition of symmetric hyperbol- 
icity for first order systems. However, there is a direct 
relationship between the Z3 systems proposed in Ref. [s^l 
and the KST formalism. This means that we can at least 
provide necessary conditions that can be helpful by re- 
ducing parameter space in the quest for a definitive an- 
swer: 

• The harmonic case is the only one in which we can 
get both pseudo-hyperbolicity and light speed as 
the only nontrivial characteristic speed. We can 
easily identify the harmonic case in the family of 
gauge conditions proposed in Ref. ^3 ; getting the 



following restrictions on their gauge parameters: 

^5 = 2, = 2(7 = 1 , £5 = -1 , EL = , A = -1 . 

(53) 

As far as the damping terms containing {ks , kl} 
are not relevant for the hyperbolicity analysis, this 
completely fixes the gauge parameters. 

• As seen in the preceding sections, the Light cones 
can be always obtained independently of the de- 
tails of the lapse and shift cones. This means that 
the results presented in Ref. [sOl for the zero shift 
case, where the non-trivial characteristic speeds 
were given explicitly in terms of the KST param- 
eters, can be applied as such. It follows that the 
light speed condition imposes the following addi- 
tional restrictions on the KST parameters: 

?7 = 4(l + x), x = 27(H-x)- (54) 

This means that ( and 7 (which amounts to the pa- 
rameter n in H43|l : n = —27) are the only remaining 
(independent) KST parameters , although ^ = — 1 
is usually required for symmetric hyperbolicity |26j | . 

There are, of course, the 10 additional coupling param- 
eters -0 introduced in Ref. [2^, so that the question can 
not be completely solved here. Nevertheless, we hope 
that the conditions we provide, which actually reduce 
by half the 22 parameter space, will pave the way to a 
definitive answer. 
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